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Abstract 

We reduce the necessary and sufficient biseparable con- 
ditions of the four qubit cluster diagonal state to con- 
cise forms. Only 4 out of the 15 parameters are proved 
to be relevant in specifying the genuine entanglement of 
the state. Using the relative entropy of entanglement as 
the entanglement measure, we analytically find the gen- 
uine entanglement of all the four qubit cluster diagonal 
states. The formulas of the genuine entanglement are of 
five kinds, for seven different parameter regions of entan- 
glement. 



Oh 
I 



> 

^' 

o 



1 Introduction 

In quantum information and quantum many-body 
physics, multipartite entanglement is still a phenomenon 
that is poorly understood [T] [5]. Experimentally, multi- 
partite entanglement has been observed in ion traps [3], 
photon polarization [3] , superconducting phase or circuit 
qubit systems [5], and nitrogen- vacancy centers in dia- 
mond [6]. More than ten qubits have been entangled 
in Greenberger-Honre-Zeilinger (GHZ) state in ion traps 
[3], and six-photon cluster state has been observed [7]- 
Usually the multipartite entangled states observed in the 
experiments are graph states. Graph states are based 
on graphs and are very useful in constructing quantum 
error-correcting codes. Due to noise and imperfections 
in preparation, the states prepared are usually mixed 
states. Typically they are the so called graph diagonal 
states. Theoretical research has been concentrated on 
the separability and biseparability of the prepared states 
[8] [9] [lO] [H] [IS [13]- Cluster states are special graph 
states. Recently, necessary and sufficient biseparable cri- 
terion for four qubit cluster diagonal states was obtained 
[Tl] . Biseparable states are multipartite quantum states 
that can be expressed as a convex sum of the projec- 
tors of product vectors and bipartite entangled vectors 
[15]. Hence a genuine multipartite entangled state is not 
biseparable. Full separable states are those that can be 
expressed as a convex sum of the projectors of product 
vectors. So full separable state set is the subset of bisepa- 
rable state set. The criterion for full separability of a four 
qubit cluster state is not known. Thus the quantification 
of the entanglement is not available when the entangle- 
ment measure involves with the full separable state set. 
Quantifying multipartite entanglement is a difficult 



problem even for a pure multipartite state. Measures 
such as entanglement cost, distillable entanglement work 
well and have clear operational meanings in bipartite sys- 
tems. However, it is not easy to extend them to multi- 
partite systems. The relative entropy of entanglement 
(REE) is a valid measure for multipartite as well as for 
bipartite systems[16]. For a given quantum state f7,REE 
is defined as E = miiip^sep S{a \\p), where Sep is the 
separable state set, S{(t\\p) — Tr{a\og2a — crlogjp) is 
the relative entropy. In multipartite system, genuine en- 
tanglement measured by relative entropy can be defined 
by minimizing the relative entropy over all biseparable 
states. Vedral et al. [16] had proposed such a definition 
of entanglement for three-partite system. In this paper, 
based on the biseparable criterion, we study the genuine 
entanglement of four qubit cluster diagonal states mea- 
sured by REE. 



2 Cluster diagonal state and its 
necessary and sufficient bisepa- 
rable criteria 

A simple graph G — {V, F) is composed of a set V of 
n vertices and a set of edges characterized by the ad- 
jacency matrix F. The n x n symmetric matrix F has 
nullified diagonal elements and F^ = 1 if vertices i and 
J are connected and Fj^- = otherwise. The neighbour- 
hood of a vertex i is denoted by Ni = {j ^ V iTij — 1 }. 
Consider a system of n qubits, and define the mutually 
commutating stabilizer operators: 
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where Xi and Zj are Pauli X and Z matrices at vertices 
i and j, respectively. The operators stabilize the graph 
state 
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such that Ki\G) = |G) . Graph state basis are 
\Ga^a2...aJ = Z"' ' ' ' Z^" \G) with a, = 0,1, each 
of them is the common eigenstate of all the operators 
Ki, with eigenvalues ±1. 

The graph of the four qubit cluster state \Cl) has three 
edges such that F^^i+i = F^+i^i — I (^i — 1,2,3) and all 



the other elements of F are 0. The stabihzer operators 



are 



= X1Z2/3/4, K2 = ZiX^Z^h, (3) 

= /l^2^3^4, = hhZsXi. (4) 

where li is the identity matrix for vertex i. The 
four qubit cluster basis states are \Claia2a3ai) with 
Ki \Cla1a2a3a4) — 1^^^010:20:304) • ^ four qubit 

cluster diagonal state a is the probability mixture of the 
cluster basis states 
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where Fai3jS > and J2a,fs,f,s=o ^^p-iS = 1- 

The original necessary and sufficient biseparable crite- 
ria are [Tl] 



definition of p2o+(5- If (O and/or are violated, then 
at least one case of ^ is violated for some /3, 7, /i, i^, since 
and (|15p are special cases of © • The same reasoning 
can be applied to prove the equivalence of inequality 
and inequalities (fT3)) (IT4)) (IT6l) (fTTl) . 

If inequahty (fT2|) is violated, then inequalities (fT5]) - (fT7)) 
are preserved due to the normalization X^qA 7 s ^ap-yS = 
So,5(P2o+<5+P4+2o+5) 1- If inequality Q is violated, 
we also have inequalities (fTSjl - dTT]) been fulfilled. We con- 
clude that if one of the inequalities P^ - (|Ti)) is violated, 
then inequalities (fT5)) - ((T7| are all preserved and vice 
versa. So we only need to consider half of the inequali- 
ties been violated. We consider the violation of inequal- 
ities (fT^ - P^ . in the following we will mainly work on 
the parameters po^ps^pi^pr. Suppose, for example, the 
maximal of -fo/370 Wil = 0, 1) be Foooo, the maximal of 
-F1/37I il3, 7 = 0, 1) be Fiooi, then po = i^oooo , P3 = -Fiooi, 

P4 = -foOlO + ^0100 + -P0IIO7 P7 = -FlOll + -FllOl + -full- 
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2Fai3^s + '^F-^^g < {Fa^^s + F^^^j 

+F^^„s + F-^^g) (7) 

for all the subscripts a,/3,7, 5, /i,i^ = 0,l, where Fai3-ys = 
{Ciap-ysl c |C4o/37(5) ■ Violation of the biseparable criteria 
means genuine entanglement. Note that the right hand 
side of inequality ([7]) is just equal to 1 according to the 
normalization of a in cluster state basis. Hence the two 
criteria can be written as 
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Denote 



P2a+S = max{i^Q/3^i-}, 

P,7 



P4,+2a+S 
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Fap-yS — P2a+S- 



(8) 
(9) 
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Then inequalities ([U and ([9]) can be further reduced to 
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We will prove that inequality ^ and inequalities (fT2|) . 
([TS]) are equivalent. For any P.^.^i.v, if ^ is violated, 
then one of and should be violated from the 



3 Entanglement measure 

The (genuine) REE [TB] defined for a tripartite entangled 
state can easily be extended to a generic multipartite 
state. The genuine entanglement of a genuine entangled 
state cr measured by the REE is 



E ~ min S{a\\p) 

p^Bisep 
= Tr(cr log2 CT - cr l0g2 £»), 



(18) 



where Bisep is the biseparable set, g is the closest bisepa- 
rable state, namely, biseparable state that minimizes the 
relative entropy. The genuine REE measures how 'far' 
is the genuine entangled state from its nearest bisepar- 
albe state, i.e, the state that is not genuine entangled. 
For a cluster diagonal state a, it is easy to show that 
g is also a cluster diagonal state following the reason- 
ing of declaration that the closest separable state of a 
Bell diagonal entangled state is a Bell diagonal state 

[IS]. Let g = J2a,l3,-y,S=0^al3-tS\Clal3'f5) {Cla/S'ysl, with 

{Aa/375} forming a probability distribution. The genuine 
REE is 

1 
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The convex property of — log function means that the 
closest biseparable state g should be at the bound- 
ary of the biseparable state set. Denote X2a+s — 
max^_^{AQ/3^5}, A4+20-1-5 = S/3,7 Ao075 — A20-1-5, then at 
least one of the equality should be reached in the bisepa- 
rable criteria of the state namely, if we replace pi with 
Xi for inequalities p^ - pTl) . then at least one of the in- 
equalities should be an equation. We use pi to specify 
the possibly genuine entangled state cr, use Xi to specify 
the closest biseparable state g. 



4 Three parameter states 

Suppose a genuine entangled state cr violate inequality 
p^ or inequality P^ or both of them, while leaving 
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Figure 1: (Color online) Edge lines or curve and gen- 
uine entanglement or biseparable regions for po — 0.3. 
Genuine entangled state Regions: Ai violates ([T^ and 
([T3)) . A2 violates ((T2]), the genuine REE can be Ea,Ea" 
or Ea'" due to the fourth parameter p4; B violates (IT^ 
and (HHD, the genuine REE is Eb] Ci violates (jl2p and 
P^ . C2 violates (|13l) . the genuine REE is Ec- Bisepara- 
ble state region: D2- Region whose biseparablity yet to 
be determined by the fourth parameter: Di. In regions 
Di,D2, both and ([T^ are preserved. For all the re- 
gions, possible violation of ([H)) will be discussed in the 
next section. 
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Figure 2: (Color online) Edge lines or curve of genuine 
entanglement regions for po = 0.6. Genuine entangled 
state Regions: A,B,C, both (fT^ and ([T^ are violated. 
The genuine REE for region A can be Ea , Ea" or Ea'" 
due to the fourth parameter p4, The genuine REE is Eb 
for B, is £;c for C. 



inequality ([T^ pending for further consideration in the 
next section. Hence a is characterized by three parame- 
ters Po,P'i,P7 concerning its genuine entanglement prop- 
erty. Without considering the fourth parameter p4 or 
possible violation of inequality ([T4|) . the genuine entan- 
glement obtained in this section is only the candidate of 
the final result. 

When inequality (IT^ is violated, namely, pa+py, > 5, 
the parameter regions can be shown as Ai, A2, i? and Ci 
in Fig.l for po < ^ or A, B, C in Fig. 2 for pq > \. When 
inequality (fT3)) is violated, namely, 2po +P3 +P7 > 1, the 
parameter regions can be shown as A2, B, Ci and C2 in 
Fig.l for Po < 5 or A, B, C in Fig. 2 for Po>\. 

We first consider the case of po < 5- For any give 
genuine entangled state in regions Ai,A2,B, Ci and C2, 
there are three classes of possible closest biseparable 
states. 

Class I is the closest biseparable state set with 



Ao + A3 = I. 



(20) 



Using Lagrange multiplier to minimize the relative en- 
tropy, we can determine 



Ao - 



Po 



2(P0+P3) 



, A3 = 



P3 



(21) 



2(1-P0"P3) 

and the genuine REE is 



2(P0+P3) 

for aU the others (22) 



Ea=^1- H2ipo+P3) 



(23) 



where _ff2 (2;) — — a; log2 a;— (1 — x) log2(l — a;) is the binary 
entropy function. 

Class II is the closest biseparable state set with 

2Ao + A3 + A7 = 1. (24) 

We obtain the solution 

Ao = ^(1-P3-P7), A3=P3, X7=P7, (25) 
(1 - P3 - P7)Fai3jS 

2(1 -po -P3 -P?) 
and the genuine REE 

Ec = (l-p3-P7)[l-ff2( 
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for all the other^G) 



Po 



1 - P3 - P7 
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(27) 



Class III is the closest biseparable state set with both 
(|^n|) and ([M)) . The Lagrange multiplier solution of min- 
imizing ([T9| with restrictions of (|20)) and ([24| and the 
normalization condition is 

Ao = ^(1-P3-P7), A3 = A7 = i(p3 +P7J28) 



A, 



(1 - P3 - P7)FaP',S 

2(1 -po -P3 -P7) 
The genuine REE is 

Eb = I - {P3 + P7)H2{- 



for aU the others(29) 
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For class I closest biseparable states, we should check if 
2Ao + A3 + Ay < 1 is preserved. Otherwise, the supposed 
closest state is not biseparable, the solution (l?T|) - (l^^ is 
no longer proper. Using the solution (j21D - (P^ . we have 



2Ao + A3 + A7 



1 



than or equal to 1 when 



P7 



1-P0-P3 



), which is less 



1 - Po - P3 

P7 < \ P3 = PAB- 

P0+P3 



(31) 



This characterizes the region A (the union Ai and A2 in 
Fig.l). The equality in (PT|) represents the border curve 
in Fig.l or Fig. 2 that divides region A and region B. 

For class II closest biseparable states, we should check 
whether Aq + A3 < ^ is preserved. If it is not, the sup- 
posed closest state will go out of the biseparable state 
set, the solution (pSjl -ipT ]) is not valid. Using (pS)) . we 
have Ao + A3 = 5(1+^3—^7), which is less than or equal 
to i only when 

P7>P3- (32) 

This characterizes region C (the union Ci and C2 in 
Fig.l). The equality in (j32l) represents the border line 
in Fig.l or Fig. 2 that divides region C and region B. 

All the states in region A violate ((5^ . thus they have 
not closest biseparable states of class II. The possible 
solutions are class I and class III with genuine entan- 
glement Ea and i^s, respectively. Since 1 — H2{x) is a 
convex function, we have Eb > Ea, the genuine REE of 
region A is Ea- 

All the states in region C violate ([51]) , so they have not 
closest biseparable states of class I. The possible solutions 
are class II and class III with genuine entanglement Ec 
and iSs, respectively. However, Ec is only a part of Eb, 
we have Eb > Ec- Hence, the genuine entanglement of 
region C is Ec- 

The states in region B (the states at the border curve 
with A and border line with C are not included) violate 
either ([?!]) or ([5^ . Their closest biseparable states can 
not be in class I or class II. They are in class III. The 
genuine REE of region B is Eb- 

The case of po ^ 5 can be similarly analyzed. In Fig. 2, 
all the states in regions A, B, C are genuinely entangled. 

The conclusion for three parameter states is: the can- 
didate genuine REE in regions A, B and C is Ea, Eb and 
Ec- Note that p3 is the maximum of Fi^^i, so 3p3 > py. 
Hence the region with < ^pj is meaningless. 

The closest biseparable states of the genuine entangled 
state in region A and C are on the lines of (PH)) and 
(hyperplanes in fact ), respectively. The closest bisepa- 
rable states of the genuine entangled state in region B 
are on the intersection point of the lines ([20)) and ([M]). 



5 Four parameter states 

The regions are classified as A,B,C,D with parameters 
Po,P3,P7 in last section. We will consider the further 
classification of each regions by the fourth parameter p4 . 
The biseparable state set is shown in Fig. 3 for Aq = 0.2. 
The three classes of closest biseparable state sets defined 



in last section are shown in Fig. 3 with surfaces I, II and 
intersection line III. There are two new classes of closest 
biseparable states appear when the fourth parameter p4 
is considered. 

Class IV (shown in Fig. 3 with surface IV) is the closest 
biseparable state set with 



An + 2A3 + A4 = 1. 



(33) 



The Lagrange multiplier solution of minimizing (jl9p with 
restrictions of is 

A3 = ^{l - po - Pa), Xo=Po, A4 = P4, (34) 

Aafl-yi — — — — for all the otherfc35) 

2(1 -po -P3 -P4) 

The genuine REE is 



Ea'" = {I - Po - Pdi^ - H2iY 



P3 



PO-Pi 



-)]. (36) 



Class V (shown in Fig. 3 with line V) is the closest 
biseparable state set with both (PU)) and ([55t . The solu- 
tion of (HH) is 

A3 = i(l-po-P4), Ao = A4 = ^(po +P4i37) 



A, 



a/37(5 



(1 - PO - P4)Fap^s 

2(1 -po-P3 - Pa) 
The genuine REE is 



for all the others(38) 
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PO 



^(1 -PO -P4)^^2( 



PO +P4 
P3 



1 - PO - P4 



)]• (39) 



The genuine entanglement regions in the four param- 
eter system with given po are determined by parameters 
, , 1-2^ , . All the genuine entanglement regions 
described in Fig.l and Fig. 2 now are three dimensional 
when p4 is considered. The regions are three dimensional 
by adding the third dimension p4 to the two dimensional 
graphs shown in Fig.l or Fig. 2. For example, the bottom 
of the three dimensional region A is the graph in Fig.l 
or Fig. 2. The roof of A is determined by the condition 



Po +P4 +P3 +P7 = 1, 



(40) 



which comes from the normalization condition 1 = 
SQ,/i,7,5 ^"/37<5 > Po +P4 +P3 + P7- The othcr border 
surfaces are P7 = 0, py = pab and Po + P3 = |. 

For region A, we have proven at last section that the 
closest biseparable states can not be in class II. Suppose 
the closest biseparable states are in class I, the solution is 
m-m- Then Ao + 2A3 + A4 ^ i(l + ^ + j^^)- 
The biseparable condition Aq + 2A3 + A4 < 1 is preserved 
only when 



1 - Po - P3 
P4 < ; Po = PA'A"- 

P0+P3 

We denote the region in A limited by (|1T|) as A'. 



(41) 



For the other part of A, suppose the closest biseparable 
states be in class IV, we have the solution Then 
'^o + ^3 = ■^{^ + Pq — Pi), which is less than or equals to 
i when 

Pi>Po- (42) 

We denote the region in A limited by ([^^ as A'". In 
region A, we have po + Ps > \, from which we can derive 
Po > PA'A"- Hence A' and A'" do not overlap. The region 
in A with po < P4 < PA' A" then is denoted as A" . Hence 
the region A is divided into three layers from bottom to 
top as A' ,A" and A'" along p4 direction. The layers are 
shown in Fig. 4 without considering parameter p^. 

For layer A" , the closest biseparable states neither be- 
long to class I nor class IV. We consider class V. The 
solution then is ([SZl)-®- It is possible that A' and A'" 
may have class V closest biseparable states. However, 
we have Ea < Ea" due to the convexity of the function 
1 — H2{x). We have Ea'" < Ea" for is a part of 

Ea'i- Thus closest biseparable states of A' and A'" can 
not be in class V. 

We should further check if the closest biseparable 
states of A are in class HI. We have proven that Ea < 
Eb in last section, the proof is valid for layer A' . We 
only need to consider A" and A'" . The closest bisepa- 
rable states for layers A" and A'" can not be class HI 
(see appendix). The genuine REE is Ea^Ea" tEa'" for 
A' , A", A'", respectively. 

For all states in region A, inequality (|12p is violated. In 
layer A', inequalities (IT5|) and (ITil) can be violated or not 
due to the location of the state. In layer A", (jl3l) may 
or may not be violated due to the location of the state. 
In layer A'" , (fT3|) is preserved since 1 > po+Pi+pz+pj 
and p4 > po- In both layers A" and A'", is violated 
as shown in Fig. 4. 

Consider the closest biseparable states of class IV, we 
have p4 > po in order to preserve Aq + A3 < \. When ([T^ 
is violated, as in regions B and C, we have 2po+P3+P7 > 

I > Po +754 +753 the last inequality comes from the 
normalization condition. Thus po > Pa- So regions B 
and C have not class IV solution. Class V is also not 
a solution for them (see appendix). The genuine REE 
keeps the same form as in the three parameter situation 
of last section for regions B and C. In regions B and C, 
inequality is violated, in regions B and Ci , inequality 
([T^ is violated. Inequality can either be violated or 
preserved in region B. In regions C and D2, the inequality 
([T4)) is preserved. Note that po +P4 +P3 +P7 < 1) we 
have Po + 2p3 + p4 < 1 + P3 — P7, thus if P3 < P7, we 
should have inequality (ITil) . Region D2 is biseparable 
even considering the fourth parameter p4. 

Violation of inequality is possible for regions Di as 
shown in Fig. 4. Region Di is divided into two layers D'^ 
and D'l, with D'( biseparable and D'^ genuine entangled. 
Consider the class VI closest biseparable states for layer 
D'l , the genuine entanglement of layer D'^ is Ea'" ■ Closest 
biseparable states in class I can not be the solution since 
D'l does not overlap with A' . Closest biseparable states in 
class V can not be the solution since Ea'" < Ea" ■ Class 

II solution requires P7 > ps, while in D[ we have P7 < P3. 




Figure 3: (Color online) Border surfaces of biseparable 
state set for Ao = 0.2. Surfaces I, II, IV are for Aq -I- A3 = 
1/2, 2Ao + A3 + Ay = 1, Aq + 2A3 + A4 = 1, respectively, 
they are the border surfaces of biseparable and genuine 
entangled states. Surfaces VI and VII are for Aq -I- A3 -I- 
A4 -f- A7 = 1 and A4 — 3Ao, they and the surface A7 = 
3A3 are the border surfaces of biseparable states and the 
unphysical region. Lines HI and V are the intersections 
of the surfaces. 

Class HI is also not the solution for D[ (see appendix). 

The situation of p4 > 3po does not exist due to the 
assuming of po = max^^^ -fo;370- 

The parameter regions and the corresponding genuine 
REE are summarized in Table 1. All the regions should 
be subjected to the normalization constrain p4-|-po-|-p3 -I- 
P7 < 1 and p4 < 3po, P7 < 3p3. 

6 Conclusions and Discussions 

We have completely solved the genuine entanglement 
problem of four qubit cluster diagonal state. For any 
probability mixture of four qubit cluster basis states, we 
first reduce the total number of parameters from 15 to 
8. Then we have proven that at most half of the nec- 
essary and sufficient biseparable criteria of are violated 
for any genuine entangled state, the number of the pa- 
rameters involved then is further reduced to 4. The four 
parameter states are classified as biseparable and genuine 
entangled. The entanglement measure for a genuine en- 
tangled state a is the relative entropy of a with respect 
to the closest biseparable state. We found that the clos- 
est biseparable state set is on the interface of biseparable 
state set and genuine entangled state set. The closest 
biseparable state set is divided into five classes. We clas- 
sify the genuine entangled states into several regions and 
find their closest biseparable state classes. Each region 
has its closest biseparable state set. The entanglement 
is given analytically for any genuine entangled four qubit 
cluster diagonal state. 

The four parameters are symmetry in some sense. For 
Po + P3 > ^! the parameter region is divided into five 




Figure 4: (Color online) Edge lines or curve of genuine 
entanglement regions for po — 0.3. Genuine entangle- 
ment regions A', A" are based on the regions Ai or A2, 
region A'" is based on A2 in Fig.l, further character- 
ized by the fourth parameter p^. The entanglement is 
Ea,Ea",Ea"', respectively. Biseparable region D'{ and 
genuine entanglement region D[ are based on region Di 
in Fig.l, further characterized by p^. The dash line is the 
border of inequality . 



The genuine entangled or biseparable regions 





Descriptions 


E 


V 


A' 


P4 < PA'A";P7 < pab; 
P3> \-Po 


Ea 


(HI* 


A" 


PA' A" <Pi< po; Pi < pab; 

P3> \-P0 


Ea" 


m 


A'" 


P4>Po; P7<pab; 
Pi,> \-Po 


Ea'" 


m 
m 


B 


P3> PT> pab; P3> -^-Pa 


Eb 


eg) (USD 

d* 


Ci 


P7 > P3; P3> i-PO 


Ec 




C2 


P7 > 1- 2po-p3; 

P3 < 5 - Po 


Ec 




D[ 


P4 > 1 - po - 2p3; P7 < Ps; 

P3<k-P0 


Ea'" 


m 


D'l 


P4 < 1 - po - 2p3; pr < pz; 

P3 < 5 " Po 


0* 


None 


D2 


1 - 2po -P3 > P7 > Pa; 

P3 < i - Po 


0* 


None 



Table 1: The regions classified by their parameters. Here 
E is the genuine entanglement. V means violation of 
inequalities ([ni)(II31)and (H!]). (dS])* or (HH)* represents 
that (Hni) or (fTi)) is violated in some part of the region 
and preserved in the other part of the region. 0* is only 
when (fT5|) - ([T7)) are preserved. 



subregions: A'" , A" , A' , B and Ci, where Ci and A'" are 
symmetric, B and A" are symmetric. A' is self-symmetric 
under the exchange 



(P0,P4) <^ (P3,P7) 



(43) 



For Po -l-p3 < ^, the parameter region is divided into four 
subregions: D[,Di,D2 and C2. The regions D'^and C2 
are symmetric with each other in the sense of the (|43l) . 
D'l and D2 are the biseparable subregions in the four pa- 
rameter system. Notice that when all of the inequalities 
p^ - P^ are preserved, we should consider the violation 
of inequalities (jl5D - PT)) . A similar analysis should be 
added for the parameters Pi,P2,P5,P6- The actual case 
is that Z?" and D2 should be further divided into truly 
biseparable subregions and genuine entangled subregions 
according to the parameters pi , p2 , ps , pe • 

The genuine relative entropy of entanglement is an- 
alytically expressed as five formulas, according to the 
subregion the genuine entangled state belongs to. The 
five formulas can be further classified as three kinds: one 
symmetric formula and two pairs. The pair formulas can 
be interchangeable under the parameter exchange (j43l) . 
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Appendix: Comparison of Eb and 

Ea" 

Define a function 

Eix) = 1 - xi/2(-) - (1 - x)H2{-^), 

X I — X 

then Eb and Ea" can be expressed as Eb — E{xb), 

Ea" = E{XA") with XB = P3+ P7, XA" = 1 - Po - P4- 

Since P3 + P7 + Po + P4 < 1, thus 

XB < XA"- 

The derivative of the function E{x) is 



dE(x) , X 

= - l0g2 



l0g2 



X-P3 



dx 1 — X ' 1 — a; — Po ' 

which leads to the solely extremal point 

P3 

X . 

Po + P3 " 

At the extramal point x*, the function E(x) reaches its 
minimum E{x*) since the second derivative aX x = x* 
is positive. The function E{x) monotonically decreases 
with X for X < X* , E{x) monotonically increases with x 
for X > X* . 

Let the genuine entangled state be at regions B, C, we 
have P7 > PAB ■ We can rewrite it as pa -|- P7 > p^^p^ ■ 



Hence x* < Xb < xa"- At the right side of x*, the func- 
tion E{x) is a monotonically increasing function, so that 

Eb < Ea"- 

In region the genuine entanglement is Eb- In region C, 
the genuine entanglement is Ec since Ec < Eb < Ea"- 
Consider the regions A", A'" and D[, we have p4 > 
PA' A"- We rewrite it as po +P4 > — , further we have 
1— Po H — P4 < p^+p^ 1 which is xa" < x*. Hence xb < 
xa" < a;*. At the left side of x*, the function E{x) is a 
monotonically decreasing function, so that 

Eb > Ea'i- 

In regions A" , the genuine entanglement is Ea"- In re- 
gions A'" and D'^ , the genuine entanglement is Ea'" since 
Ea'" < Ea" < Eb. 
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